Vortices in small superconducting disks 
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We study the Ginzburg-Landau equations in order to de- 
scribe a two-dimensional superconductor in a bounded do- 
main. Using the properties of a particular integrability point 
(k — l/\/2) of these nonlinear equations which allows vor- 
tex solutions, we obtain a closed expression for the energy of 
the superconductor. The presence of the boundary provides 
a selection mechanism for the number of vortices. 

A perturbation analysis around k = l/\/2 enables us to 
include the effects of the vortex interactions and to describe 
quantitatively the magnetization curves recently measured on 
small superconducting diskil. We also calculate the optimal 
vortex configuration and obtain an expression for the confin- 
ing potential away from the London limit. 
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The dimensionless Ginzburg-Landau energy functional 
of a superconductor depends on only one parameteiH, 

the ratio k between the London penetration depth A and 

the coherence length ^ 

The order parameter is dimensionless as well as the 
magnetic field B measured in units of ^^p- , with (t>a — 
Lengths are measured in units of A\/2- The expression 
(|l|) assumes that both the order parameter and the vector 
potential have a slow spatial variation. The integral is 
over the volume — TiE?d of a thin disk of radius R and 
thickness d. 

Outside the superconducting sample, the order param- 
eter vanishes and the magnetic field is solution of the 
Maxwell equation. The boundary condition between a 
superconductor and an insulator is (V — iA)ip\fi = 
where n is the unit vector normal to the surface of the 
disk. The presence of a boundary precludes a complete 
analytical solution of the 3d Ginzburg-Landau equations 
for a thin disk. We arc then led to make some simplifying 
assumptions, based upon numerical resultaj. 

The thickness d of the sample considered in the 
experimental fulfills d <C ^ and d < \. If the curvature of 
the magnetic flux hues, given by i?/Ag (where Ae(d, i?, A) 
stands for the effective screening length) , is smaller than 
1/ Ae , i e. if 7? <gC Ae , then both ip and the vector potential 
A can be considered to be constant accross the thickness 
and the disk is effectively two-dimensional. The expres- 
sion of the effective screening length XeidjJi,X) is not 
known, except for the case R ^ oo wherdS'El Ag — A^/d. 
In the London limit {i.e. k ^ oo), |Such a system has 
been described using Pearl's solutionB. Finally, since ip 



and A are constant over the thickness, the covariant Neu- 
mann boundary condition, stated above, is automatically 
satisfied on the upper and lower surface of the disk. 

The Ginzburg-Landau equations arc nonlinear, second 
order differential equations whose solutions are usually 
unknown. However, for the special value k — A= known 

, the equations for ip and A reduce 
to first order differential equations and the minimal free 
energy can be calculated exactly for an infinite plane. 
This relies on the identity true for two dimensional sys- 
tems |(V - iA)ip\'^ = + y X J + SIV'P where j 
is the current density and the operator V is defined as 
'D = dx + idy — i{Ax + iAy). At the dual point, the 
expression (|l|) for is rewritten using this identity as 
follows 
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+ + {j + A).dl (2) 



where the last integral over the boundary of the sys- 
tem results from Stokes theorem. 

For an infinite plane, we impose that the system is su- 
perconducting at large distance, i.e. — > 1 and J — > 
at infinity so that the boundary term in coincides 
with the London fluxoid. It is quantized and equal to 
§dn ^X-dl — 27rn, where x is the phase of the order pa- 
rameter. The integer n is the winding number of the 
order parameter 7p arid as such is a topological charac- 
teristic of the systemB. The extremal values of JF, are 
J- — 27m, and are obtained when the bulk integral in (|^) 
vanishes identically giving rise to two first order differ- 
ential equations. These two equations can be decoupled 



inear equation which 
. However, for the in- 



to give for a second order noi. 
admits families of vortex solutionsi 
finite plane, there is no mechanism to select the value of 
n, which only plays the role of a classifying parameter. 

The extension of these results to finite size systems, 
namely the existence and stability of vortex solutions and 
their behaviour as a function of the applied field received 
a partial numerical answer. jMumerical simulations of the 
Ginzburg-Landau equational3 show the existence of sta- 
tionary vortex solutions whose number depends on the 
applied magnetic field. Moreover, these simulations indi- 
cate that the physical picture derived for k = remains 
qualitatively valid for quite a large range of values of k, 
with a small corresponding change of free energyllll . 

We consider finite size systems at the dual point i.e. 
for = There, the edge currents screen the external 
magnetic field therefore producing a magnetic moment 
opposite to the direction of the field, whereas vortices in 
the bulk of the system produce a magnetic moment along 
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the direction of the applied field. Assuming cylindrical 
symmetry, the current density j has only an azimuthal 
component, with opposite signs in the bulk and on the 
edge of the-system. Thus, there exists a circle F on which 
Jvanishesll^. This allows us to separate the domain into 
two concentric subdomains fJi and ^2 whose boundary 
is the circle F. Therefore one can extend to the subdo- 
main the results obtained for the infinite case. The 
existence of vorticesjn a finite domain such as l^i was 
checked numericallyo. It was shown that vanishes as 
a power law at the center of the disk, hence there is a 
(multi-)vortex in the center whose multiplicity is deter- 
mined by the exponent of the power law. The magnetic 
flux <I>(rii) — n\s, quantized and the free energy in fii is 
T(Slx) = 27rn. 

The contribution of T{VL2) to the free energy can be 
written, using the phase and the modulus of the order 
parameter V', as 



(V|V'|)' + |V'(Vx-A)P + 



+ |^|2)2 



(3) 



We know, from the London equation, that both the mag- 
netic field and the vector potential decrease rapidly away 
from the boundary of the system over a distance of 
order A-\/2. Over the same distance, at the dual point, 
saturates to unity. One can thus estimate the integral (J^) 
using a saddle-point method. We assume cylindrical sym- 
metry, and we neglect the term (VIt/iI)^ on the boundary 
because of the boundary conditions, so that the relation 
is now given by an integral over the boundary of the 
system. To go further, we need to implement boundary 
conditions for the magnetic field B{K) and the vector 
potential A(R). The choice B[K) — Be, where Be is the 
external imposed field, corresponds to the geometry of 
an infinitely long cylinder, where the flux lines are not 
distorted outside the system. A more suitable choice for 
a flat thin disk is provided by demanding (p = (jje- This 
boundary condition implies that the vector potential is 
identified by continuity to its external applied value Ae- 
It should be noticed that the magnetic field B has then 
a non monotonous variation: it is low in the bulk, larger 
than Be near the edge of the system, because of the dis- 
tortion of fiux lines, and evpatually equal to its applied 
value far outside the systemlla. 

Finally, the minimization of the free energy with re- 
spect to IV'I gives 1 — l^/ip = |Vx — A\^, such that, per- 
forming the integral over the boundary of the system, we 
obtain 

= - C^ef - \{^f{n - Pef (4) 

We have neglected the contribution of the B^ term, which 
is smaller by a factor of the order {X/Kf . 

The thermodynamic Gibbs potential Q of the system 
is then 



■T^-Qin, 4>e) = n + a{n - (j)ef - — (n - (f)e)'^ ~ a^0e^ 

ZTT Z 

(5) 

where we have defined a = The relation (^) con- 

sists in a set of quartic functions indexed by the integer 
n. The minimum of the Gibbs potential is the envelop 
curve defined by the equation §§10^ — 0, i.e. the sys- 
tem chooses its winding number n in order to minimize 
Q. This provides a relation between the number n of 
vortices in the system and the applied magnetic field 0e. 

We consider the limit of large enough such that 
the quartic term is negligible. The Gibbs potential then 
reduces to a set of parabolas. The vortex number n is 
then given by the integer part 
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while the magnetization AI = is given by 

- M = 2a{<j)e ~n)~ 2a^<j)e 



(6) 
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For tpe smaller that , we have n ~ and (— M) in- 
creases linearly with the external fiux. This corresponds 
to the London regime. The field Hi at which the first vor- 
tex enters the disk corresponds to G{n = 0) = G{n = 1), 
i.e. to 



Hi 



27rV2i?A 27ri?2 



(8) 



The subsequent vortices enter one by one for each cross- 
ing Q{n -|- 1) = G(n); this happens periodically in the 
applied field, with a period equal to AH = and a 

discontinuity of the magnetization AM = . 

There is a qualitative similarity between the results 
we derived using the properties of the dual point and 
those obtained fropi. a linearised version of the Ginzburg- 
Landau functional. However, the two approaches differ 
in their quantitative predictions due to the importance 
of the nonlinear term. 

Within the previous approximations, the expression 
captures the main features observed experimentally 
i.e. the behaviour of the magnetization at low fields 
(before the first discontinuity), the periodicity and the 
linear behaviour between the successive jumps. From 
the experimental parametersEI namely R — 1.2/im and 
A(r) = 84nm at T = OAK, we compute from our ex- 
pressions Hi — 25G and AH A.6G. These values agree 
with the experimental resultsQ to within a few percent. 
We emphasize that Hi scales like j^ , whereas AiJ^cales 
like in accordance with the experimental datald. We 
calculate the ratio of the magnetization jumps to the 
maximum value of M to be 0.20 as compared to a mea- 
sured value of 0.22. The total number of jumps scales 
like R^ and the upper critical field is independent of R 
in our theory in agreement with the experimental data. 
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At the duahty point k = l/\/2, the contribution of 
the vortices to the free energy is topological and does 
not depend neither on the precise shape of the vortices 
nor on the form of their interaction. This property does 
not hold for other values of k. For a 2d film with an 
infinitesimal current sheeit, the vortex configuration has 
been computed by Pearlcl, using the London equation, 
and differs qualitatively from the present model. Indeed, 
away from the dual point, both the shape of the vor- 
tices and their interaction modify the free energy and 
the magnetization. For instance, in the London limit 
(k — !■ cxd), radially symmetric solutions become unstable 
and different geometrical configurations of the vortices 
have different energies. This results from two contribu- 
tions to the energy, arising from the interaction between 
the vortices themselves and between the vortices and the 
edge currents. We have performed a perturbative analy- 
sis around the dual point, and we obtained that the bulk 
free energy J-'{Qi) is given by 

^T{ni) = n(l + 1-(kV2 - 1)) + I3{kV2- 1) VW(r„). 

(9) 

The part which is linear in n is independent of the po- 
sition of the vortices. It has been evaluated using a 
variational ansatal2l. The two-body interaction potential 
near the dual point is well approximated by a function 
U{rij) = U{r) where U{0) = 1, U{oo) = with /3 f In 
particular, for a configuration where all the vortices are 
close to the center of the disk, the bulk free energy is 

^H^i) = nil + ""^^ + ln\nV2 - 1) (10) 

Thus, away from the dual point, the linear term in the 
bulk free energy is not topological anymore and is modi- 
fied by the interaction. The attractive or repulsive char- 
acter of the interaction between vortices depends on the 
sign of (kv^ — 1). 

To obtain the edge contribution to the energy, we con- 
sider first the case of a single vortex placed at a dis- 
tance X (in units of AV2) from the center of the digit- 
Then, the phase of the order parameter is given bytZI 
tanx = cose- ' ^^*3re y — ax — xjR (whereas one has 
X = nQ for the symetrical case). Starting from the ex- 
pression (^), we obtain 

1 a? 
—T{^l2,v) = a{n - 4>ef ~ ~ + 2/' 

(11) 

where the function f{a,y,(f)e) can explicitely be calcu- 
lated and represents a vortex confining potential inside a 
finite superconductor for k ~ 1/ \/2. This corresponds to 
the well-known Bean-Livingston confining energy barrier 
in a 3d superconductor which has been obtained in the 



extreme type II limittSI using the London equation. It is 
important to emphasize that around the dual point, vor- 
tices are not point-like and therefore the usual expression 
of the Bean-Livingston energy barrier does not hold. 

The possible equilibrium configurations of vortices re- 
sult from the competition between the bulk and edge 
contributions to the free energy derived above. It is ei- 
ther a giant vortex at the center of the disk, a situation 
which preserves the cylindrical symmetry, or a polygonal 
pattern of small vortices. In order to evaluate the energy 
of these configurations, we generalize the relation ([l^) to 
the case of a polygonal configuration of vortices placed 
at a distance x from the center of the disk. The resulting 
energytZI is n-times the barrier contribution f{a,y,(f>e) 
obtained in ([ll| ) provided the following substitutions are 
made: a —^ na, y y"' and 4>e — *■ 

In conclusion, we have investigated the question of 
the existence and stability of vortices in small two- 
dimensional bounded superconducting systems. We have 
showr£3 that starting from the exact solution of the 
Ginzburg-Landau equations for an infinite plane and for 
the special value k = it is possible to derive an 

analytical expression for the free energy in a bounded 
system. The resulting expression provides a satisfactory 
quantitative description of the magnetization measured 
on small superconducting aluminium disks in the low 
magnetic field regime. For larger fields, we cannot ne- 
glect anymore the interaction effects due tQ_the vortices 
and the edge currents. Perturbation theoryO around the 
value K = l/\/2, has allowed us to derive an expres- 
sion for both the confining potential barrier of the vor- 
tices and the strength of the interaction between vortices. 
This provides a more refined description of the measured 
magnetization at larger fields. 
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